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We investigate the spectral function of light quarks in infinite quark matter using a simple, albeit 
self-consistent model. The interactions between the quarks are described by the SU(2) Nambu- 
Jona-Lasinio model. Currently mean-field effects are neglected and all calculations are performed 
in the chirally restored phase at zero temperature. Relations between correlation functions and 
collision rates are used to calculate the spectral function in an iterative process. 
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CN ■ I. INTRODUCTION 

c3 . 

' It is well known that short-range correlations have influence on the properties of nuclear matter and finite nuclei. The 
spectral functions observed in A{e, e'p)X 1] and A(e, e'pp)X 2] experiments are much wider spread than explained by 
mean-field dynamics. Due to short-range correlations a substantial amount of high-momentum processes is contained 
in the spectral function. 

I There have been many theoretical approaches, mainly based on nuclear many-body theory, trying to understand 
I> ■ the short-range correlations in nuclei. The results of these calculations have converged in the last years 

ill In 

, particular, the significant population of high-momentum states in the nucleon momentum distribution which is an 
' overall measure for short-range correlations is described rather well. 

Our present study is motivated by the work of Lehr et al. ^ . They have calculated the nucleon spectral function in 
I nuclear matter in a simple model based on transport theory. Direct relations between collision rates and correlation 
functions form a self-consistency problem and were used to determine the spectral function iteratively. In their 
> calculations the drastic assumption of an averaged, constant scattering amplitude was made. Hence the results are 
dominated by the properties of the available phase space. 

It has turned out that the results of Lehr et al. agree surprisingly well with experimental data and sophisticated 
_L( ' "state-of-the-art" calculations from many-body theory Q. Spectral functions, momentum distributions, occupation 
O , probabilities, and response functions are all in close agreement with other calculations. It is in particular striking that 
the only parameter in that model, the coupling strength, is enough to describe size and slope of the high-momentum 
tail in the nucleon momentum distribution in nuclear matter. 

This leads to two conclusions. First, the properties of the nucleon spectral function in nuclear matter are determined 
by phase space effects and an average strength of the short-range correlations. Second, the detailed structure of the 
interaction seems to be unimportant as long as the calculations are made in a self-consistent framework. 

Due to the success of this model, we take up the concept and make the assumption that similar conditions apply 
for the short-range correlations in quark matter. Only a few changes to the existing nucleon model are necessary to 
get a working model for quark matter. First of all, the low current masses of up and down quarks make it necessary 
to perform the calculations relativistically. This implies some complications since the relativistic, fermionic spectral 
function is a 4 x 4 matrix in spinor space. Using symmetry arguments it can be shown that a reduction to three scalar 
functions is possible. This makes all expressions look more complicated and increases the numerical efforts. However, 
it does not change the fundamental concept of the model in any way. 

In addition a sensible model for the quark interactions has to be used. Our basic assumption is that the details of 
the interaction should be relatively unimportant as long as the overall strength is correct and the relevant symmetries 
are respected. Hence we use the well known Nambu-Jona-Lasinio model (NJL model) in its SU(2) version in our 
calculations. It has the same symmetries as QCD and describes an effective pointlike interaction with a constant 
coupling strength. Usually the NJL model and its extensions are used in a mean-field approximation. In the spirit of 
our works on nucleons, we present here an approach that goes beyond the mean-field approximation in a self-consistent 
way. 
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In the present work, we use the standard NJL model and do not consider an additional attractive interaction in the 
quark-quark channel as it appears in extended NJL models that are used in recent works on color superconductivity 
The aspect of color superconductivity - albeit interesting - is beyond the scope of the present work. A more 
detailed discussion can be found in sec. ^of this paper. 

The calculations we present in the following are restricted to the chirally restored phase where the quarks are 
massless. The relativistic structure of the spectral function gets simpler in this phase, thus simplifying the analytic 
expressions and reducing the time needed for the numerical calculations significantly. We will come back to that point 
below. 

In sec. mthe theoretical aspects of the model are summarized. It is shown how expressions for the self-energy and 
the spectral function are constructed that are directly related to each other. The resulting self-consistency problem 
can be solved in an iterative approach. The technical details and the results of our numerical calculations are presented 
in sec. IIIII We show the results for the spectral function, its width and the momentum distribution of the quarks. 
The role of short-range correlations is investigated. Section llVI summarizes our results and gives an outlook. 



II. THE MODEL 



A. Basic relations 



In Ref. the model we adopt is described for nonrclativistic nuclcons. Many of the details concerning the 
definition and meaning of the basic quantities can also be found in Refs. |^ 0|. Some aspects of working with 
relativistic nucleons are discussed in Ref. The fundamental elements of our model are the one-particle Green's 

functions without ordering, and . For relativistic fermions they are defined as 



-ig<p{l,l') = (^MnV-aCl)), 



(1) 



where the V' are field operators in the Heisenberg picture. The indices denote the spin degrees of freedom; isospin and 
color indices have been supressed here. 1 and 1' denote two space-time coordinates. Since we describe fermions with 
our model the field operators are spinors and the Green's functions are 4x4 matrices in spinor space. 
At this point we also introduce the retarded Green's function g'"^^ . It is closely related to g^. 



7-^(1, 1') = e(ti - tiO g>^(l, 1') - .g<^(l, 1') 



(2) 



The Fourier-transformed Green's functions are used to define the spectral function A{p). Obviously A{p) takes over 
the matrix structure of the Green's functions 



= i[g>0{p) - g^p{p)] = -2luvg%{p), (3) 

where "Im" is given by Imi^" = -^{F — 7o-F^7o), and correspondingly Rei^ = ^{F + 7o-F^7o)- 

In addition to eq. lO the Green's functions are related to the spectral function via the energy-momentum space 
distribution function. In thermal equilibrium one has 

-^dapiP) = Aal3{p)nF{Po), (4) 
'^dapiP) = •^a0{p)[^ - npiPo)], (5) 

with the thermal Fermi distribution npipa)- At zero temperature the Fermi distribution becomes a step function, 
npiPo) = 6(t^F -Po)- 

The single particle self-energy of an interacting quantum system can be split into a mean- field part S'"^ and the 
collisional self-energies E^. Like the spectral function the self-energy has a matrix structure in spinor space: 

1') = 1') + 0(ii - ii')S>^(l, 1') + e(ii' - ii)S<^(l, !')■ (6) 

The mean-field self-energy is time local, S™^ ^ 6{ti — ti'). It describes the effects corresponding to the motion of 
noninteracting particles in a mean-field potential and is responsible for the dynamical generation of the constituent 
quark masses. The collisional part of the self-energy contains the effects that arise from decays and particle collisions 
in the medium. 
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As for the Green's functions a retarded self-energy is introduced 

E-^(l, 1') = E-f (1, 1') + e(ti - ty) [e>^(1, 1') - S<^(1, 1')] . (7) 

In momentum space the width of the spectral function is given by the imaginary part of this retarded self-energy. 
One has 

T^f,{p) = -2ImS-*(p) = z[I]>^(p) - S<^(p)]. (8) 

E™^ does not contribute to the imaginary part of S''''* {p) since it is time local. Hence the width of the spectral function 
is completely given by the coUisional self-energies S^. In Q it has been shown that iE>(p) and ~iT,^{p) are identical 
to the total collision rates for scattering into (gain rate) and out of (loss rate) the configuration {po,p), respectively. 
From analyticity it follows that the real part of the self-energy is related to F by a dispersion relation 

ReE-* (p) = E-f {p)+V ^^^^ . (9) 

J 2n Po- Po 

If the width is known over the full energy range the real part of E''°' can be calculated dispersively. 



B. Interactions in the NJL Model 



In the following, we use the standard SU(2) version of the NJL model It has been designed to resemble the 

symmetries of QCD. In particular, it allows the dynamical generation of fermion masses via the mechanism of chiral 
symmetry breaking. The Lagrangian of this effective interaction is given by 

/:n,7L = ^i^^ + G[(V;V')' + (^^TsrV-)'] , (10) 

where G is the coupling strength, independent of energy and momentum, and the are the isospin Pauli matrices. 
The small current quark masses are neglected here. 

The pointlike interaction renders this model nonrenormalizable. In our calculations we deal with this problem by 
introducing a cutoff for the three-momenta. This cutoff and the coupling constant are the free parameters of the NJL 
model. Usually their values are chosen so that the model reproduces the correct values for the quark condensate and 
the pion decay constant in vacuum. 

Before we continue with the discussion of our present model we want to comment on its relation to other present 
investigations. Recent works on color superconductivity 0, IT^ IT^ consider an extended NJL model. In the case of 
two flavors the interaction is given by 

Ant - Gi[(V;^)2 + i^lJ^f^P)^] + G2(l^l75T2A2^^)(V'^«75r2A2^), (11) 

where C is the charge conjugation operator and A2 is a Gell-Mann matrix in color space. Gi and G2 are coupling 
constants. In the mean-field approximation, the new term on the right of eq. (|ll|l introduces an attractive interaction 
in the quark-quark channel that can lead to a pairing of quarks at the Fermi surface. This pairing breaks the color 
symmetry spontaneously and gives rise to color superconductivity. We note in the passing that the use of such an 
extended Lagrangian is not free from ambiguities 0, Q| due to the freedom of Fierz rearrangements. 

As we will discuss below, we go beyond the mean-field approximation in our approach. We have decided as a 
first step to construct a consistent model based on the well established Lagrangian (|10|l . An additional attractive 
quark-quark channel is not considered at this stage. Color superconductivity is definitely an interesting phenomenon 
and can be incorporated into the model in the future. 



C. Structure of the relativistic expressions 



As shown above the spectral function of relativistic fermions has a matrix structure. The most general form of this 
structure is found by writing the spectral function in terms of the 16 linear independent products of the 7 matrices 
(Clifford algebra) and demanding the invariance under certain symmetries fisf . We consider here parity and time- 
reversal symmetry. A strongly interacting system of infinite size and in thermal equilibrium is invariant under these 
transformations. In the rest frame of the medium this leads to the following form for the spectral function: 



■^{p) = Ps(po,P^) + Pq{po,p'^) ■ 1° + Pv{po,P^)p- 7, 



(12) 



4 



where p is a unit vector in the momentum direction. As a first step we restrict our calculations in the present work to 
the chirally restored phase. Since the "scalar" term is not chirally invariant it has to be zero; the spectral function 
is then completely determined by the two scalar functions po and p^. Note that the spectral function for free quarks 
with po > is given by ^ = 7r|p|~^ ^S{pa — According to our choice of signs in eq. H12I) this means that po > 
and pv < in this case. 

The structure of the self-energies and the width F are identical to that of the spectral function. They too consist 
of three scalar functions, 

= J:f{po,P^) + 5]|bo,p') • 7° + S|(Po,p'¥- 7, (13) 

Tip) = r,(po,p') + ro(po,p')-7° + rv(po,p")p'-7, (i4) 

with = Fg = in the chirally symmetric phase. 

D. Self-consistent expressions 

The mean-field part S™^ of the self-energy is given by the tadpole diagrams. Currently chiral symmetry is restored 
"by hand" in our model. This means that the system under consideration is forced into the chirally restored phase by 
setting the constituent quark masses, i.e., the Hartree term of S™^ to zero. The Fock term contains a part that needs 
not to be zero but this part can be absorbed in a redefinition of the chemical potential (llj . Note that one can do 
this for any temperature and chemical potential since the gap equation for the quark masses Q is always solved 
for TO = 0. However, thermodynamically this might not be the favored phase when there exists also a finite solution 
for the gap equation. 

The collisional self-energies can be calculated diagrammatically. We restrict ourselves to the lowest order, given 
by the direct and the exchange Born diagram (see Fig. In this approximation only two particle correlations are 
included in the model. The corresponding processes at finite chemical potential (cf. Ref. j^) are shown in Fig. |21 

In perturbation theory the particle and antiparticle lines in the diagrams of Fig. ^are interpreted as free propagators. 
Our model, however, takes into account the in-medium character of the intermediate states. The lines are interpreted 
as full relativistic Green's functions which again depend on the self-energies via the spectral function (eqs. Q and 
This way we end up with self-consistent expressions for the self-energies. In the language of Feynman diagrams the 
use of full Green's functions means that we sum over a whole class of diagrams. Note that no ambiguities concerning 
Fierz rearrangements appear here. All interactions at the two loop level^ caused by the interaction term in eq. H1(J|) 
are included in our approach. In particular, we find contributions from quark-quark as well as from quark-antiquark 
scattering (cf. Fig. 

At this point it might be illuminating to discuss which Feynman diagrams are resummed by our approach and 
which ones are not. We use full propagators for the calculation of the sunset diagram depicted in Fig. ^ In this way 
we include diagrams like the ones shown in Fig. O We do not sum up the series shown in Fig. 21 (a) (we only have 
the first term). Series of that type would correspond to the coupling of dynamically generated pious or sigmas to the 
quarks [lTllT3 |. Clearly such series are relevant in the chirally broken phase. We also do not sum up the series shown 
in Fig. 01(b), which corresponds to coupling to diquarks. This is expected to be relevant if the phenomenon of color 
superconductivity is tackled. 

Working out the diagrams and replacing the Green's functions by spectral functions according to eqs. and 



^ Note that using full propagators instead of free ones implies that only two-particle irreducible diagrams have to be considered. 
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(O one gets for the two components of iT,^{p) in the chirally restored phase 

G2 



(27r) 



14 



dkodqo j dv^dv^ j d cos duvl^ — np{k())]uv 
<A 



Poik)paiq) - Pv(g)Pv(fc) 



\u + v\\u — v\ 



"F(go)[l - npiro)] 



+v^+2\p\v 



2\p\v yp.2+„2_2|p-| 

Po{k)po{q) - py{q)pv{k) 



u + v\\u ~ V 

2 _ ^„ „2 



[1 - nF{qo)]nF{ro) 



2|P|U Jp^+u^-2\p\i 



fc=l(a-ir) 

'■0=90 + '=0 -PO 



(15) 



and 



G2 



(27r) 



14 



dkodqo j du dv j dcos-duvi^ — 'n'F{ko)]uv 
<A 



(16) 



PQ{k)po{q) - Pviq)Pvik) 



\u + v\\u ~ v\ 



"-f('7o)[1 - npiro)] 



■.,-1 , dr^p^{rQ,r^)- 



''0 = P0 + 90-'=0 



Pa{k)pa{q) - Pv{q)Pv{k) 



? 2 



|m + w||it — 



[1 - nF{qo)]nF{ro) 



dr py{ro,r )- 



1 I — u 



r \ 



J ,-=i(fl+i?) 

'■0 = '!0 + '=0-P0 



Here u — q + k and v — q — k are the center of mass and relative momenta in the scattering processes of Fig. |21 "ffuv 
is the angle between the two vectors u and v. The subscript " < A" at the integrals indicates that a three-momentum 
cutoff is applied. This notation requires some comments: The integrals in eqs. H15() and (|16() are integrals for u and 
V. In practice, however, the cutoff is applied to k and q'm the calculations. The vectors k and q are the momenta of 
the particles and antiparticles in Fig. ^ thus they must be regularized when working out the diagrams. 

The expressions for — zEq (p) and —iY,^{p) are not presented here. They are easily found from eqs. p5|) and 
by replacing all functions np by {1 — uf) and vice versa. 

Eqs. p5|) and p6|l relate the coUisional self-energies (collision rates) to the spectral function. Note that the 
zeroth and the vector component of eq. H12|l are mixed. Finally, the components of the spectral function have to be 
expressed in terms of the self-energies via the width F. The expressions are found by inserting the explicit form 
of the relativistic in-medium Green's function. 
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(17) 



into the definition of the spectral function . Since the chirally restored phase is considered the quarks are massless 
and no mass term appears in the denominator. 

When working out the imaginary part of 17''°' one also needs the real part of the self-energy E''''*. In principle it can 
be calculated from F by use of the dispersion relation . The studies of Lehr et al. have shown that the influence 
of ReE''''* on the nucleon spectral function and the nucleon momentum distribution in nuclear matter is small. In 
their first calculations they have replaced ReE'^'''(p) by a constant value independent of energy and momentum. It 
was found that this violation of analyticity leads only to minor effects in the results. 

In our present calculations, we are not able to use the dispersion relation Q due to technical reasons. We work on 
a finite grid in energy and momentum space, the boundaries of which have been chosen to include all significant parts 
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of the spectral function. The width Fq, however, extends much further into the high energy regions than the spectral 
function. As we do not know Fg outside the grid, the pq integration of the dispersion integral cannot be performed^. 

In principle, one could find ReE''^* in a different way. Using the dispersion relation for the retarded propagator, 
Reg'"* can be calculated from the spectral function. Then it is possible to deduce ReS™* from eq. p7ll . We have 
checked this and it turned out that the real part of the self-energy is very small compared to po and \p\, i.e., in ^— ReE'^"* 
the dominant part is ^. Since we know from the nucleons that the real part of S''''* is relatively unimportant Rel]'^°* 
is set to zero in the following. 

The following expressions for po and pv are found from eqs. and (|17|) . when ReS'^^* is set to zero 

, .2, _ jpl+p'Wo + i(rg - r^)ro + 2po|piF. 
^ - b^-P--ir§ + iF?p + [poro-F.|pT' 

^ ~ [Pl-P'-\TI + \VW + [P.T,-TMY- ^^'^ 

The two relations look very similar. Up to some signs only Fq and Fv are exchanged in the numerator. Near the 
on-shell point po and pv differ only by a minus sign'^ at positive po if terms ^ 0{T^) are neglected 

2^ 2pgro + 2po|po|F^ 2N _ 2p^F, + 2po|po|Fo 
Po{Po,Po) ~ — > Pv[Po,Po) ~ — • (20) 

Eqs. 1)1511 and l|16() . (|18|l and (|19|l . and eq. lO form a set of equations describing a self-consistency problem. A direct 
solution of this problem is not possible but the equations can be used to calculate the spectral function iteratively. 

III. NUMERICS AND RESULTS 
A. Details of the calculation 



All calculations were performed at zero temperature in the chirally restored phase. We have used two different 
quark densities. Motivated from the investigations for nucleons in nuclear matter, we chose the quark matter in the 
first case such that it is comparable to regular nuclear matter. As every nucleon consists of three valence quarks the 
quark density was set to = 3 • pnm — 3 • 0.17 fm^^. For the number of flavors we use A^f = 2 (up and down). This 
yields a Fermi energy and a Fermi momentum of top — Pf — 0.268 GeV for massless quarks. In reality quark matter 
with that chemical potential would not be in the chirally restored phase. Therefore, we have also worked with a three 
times higher density, — 1.53 fm^"^, corresponding to a Fermi energy of ujf = 0.387 GeV. This case seems to be 
more realistic for chirally symmetric quark matter, the chemical potential is well beyond the chiral phase transition 
in the NJL model [Tllll3|. Note that all states at energies up to are filled. This includes the states at negative 
energies. Thus there are no holes in the Dirac sea which could be identified with populated antiquark states. 

The three-momentum cutoff A = 653 MeV and the coupling constant G = 2.14 • A^^ of the NJL model were chosen 
so that the model gives the known values for the quark condensate, {uu) = (dd) = — (250MeV)'^, and the pion coupling 
constant, = 93 MeV, in vacuum 'Tl'|. In the case of the lower density we performed also calculations with two 
times and four times larger coupling strengths to investigate the influence on the spectral function. 

The calculations were carried out on an energy-momentum grid with boundaries —653 MeV < po ^ 653 MeV and 
\p\ < 653 MeV (in correspondence to the NJL cutoff) and a mesh size of 6.53 MeV in both directions. Due to the 
cutoff a grid of that size is sufhcient to include all significant parts of the spectral function. The states with negative 
energies po are interpreted as anti-quark states with positive energies \po\ in the discussion of the results. 

To initialize the calculations a constant width, Fq = 1 MeV and Fv = 0, was used. In principle it would be 
interesting to start from a quasiparticle approximation with a width much smaller than the final result. This would 
allow to observe the redistribution of strength away from the peaks during the iterations. However, this is numerically 
not feasible. So the initialization values are not too far away from the final results and only two iterations were 
necessary to reach self-consistency. 



^ Keep in mind that the NJL cutoff is a three-momentum cutoff. While all functions are zero at momenta above the cutoff there is no 
direct consequence for high energies. In fact the cutoff does affect the high energy behavior of Fq (cf. eqs. 1151 and 1161 ) but only in 
the form of a contiuous decline (we will come back to that point at the discussion of the results). 

^ Concerning the sign of pv recall our remark after eq. 1121 . 
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B. Results 

In this section, we present the resuhs of our iterative calculations. Before the spectral function is discussed we have 
a look at the width T. This is helpful to understand the resulting spectral functions and it allows a clearer observation 
of the iterative process than the spectral function itself. 

We show the results for po and Fq but not for and Fv. Most of the interesting structure of the spectral function 
is located around the on-shell peaks. It has been discussed that po and p^ differ at most by a sign near those peaks 
(eq. H2()|l V Thus it is not surprising that both look very similar and it is sufficient to show only one of them. The 
reason to show po is that the density of states (particles minus antiparticles) is given by po alone, 

Tr{7o^(p)} - 4 • NiN,poipo,P^). (21) 

There is no simple relation between Fq and Fv justifying to show only one. But our results show that Fv is much 
smaller than Fq (this has also been observed in calculations for relativistic nucleons Hence Fv seems rather 

uninteresting since its influence on the spectral function is minimal. 

Figure El shows cuts of the width Fq at several momenta for uip = 0.268 GeV. The dotted line is the result after 
the first iteration while the solid line displays the final, self-consistent result. Physically the most interesting area lies 
in the energy range < po < (^f since that is the location of the populated quark states. All states above the Fermi 
energy as well as the antiquark states at negative po are unoccupied (no holes in the Dirac sea). Fq takes on low 
values of 0.1 — 1 MeV in this zone, being close to the initialization value. 

The two most important features of Fq are clearly visible. First, the width is zero along the Fermi energy ujp, 
the quarks along this line are quasiparticles. Due to Pauli blocking and energy conservation it is not possible to 
scatter into or out of this configurations at zero temperature. Second, the width grows explosively at high \po\. The 
reason for this is the pointlike interaction of the NJL model. The self-energies essentially sum up the phase 
space available for scattering processes. The opening of this phase space at large \po\ can be directly read off from 
Fg. The three-momentum cutoff applied in the calculations stops this inflational behavior at higher energies. These 
energies, however, are outside of our grid. A simple (on-shell) estimate shows that the maximum width is reached at 
A < \pci\ < 2A and that the width must be zero for \po\ > 3 A . 

Another noticeable feature is the local minimum in the antiparticle sector, most visible at low momenta. This is a 
phase space effect that can be easily explained by the energy and momentum dependence of the particle correlations. 
In the bottom row of Fig. |21 all processes contributing to the width at negative po are shown. Phase space opens 
for the processes pk — > qr and p — > kqr at po = and po = —uip, respectively, and grows with decreasing pq. On 
the other hand the process pkq ^ r is possible for po > —\p\ only. Its phase space grows with increasing pq and is 
maximal at p^ = 0. The overlap of these three processes leads to the observed minimum that is situated approximatly 
at Po — — \p\ for low momenta. 

In Fig. El the width Fq is shown for ujp = 0.268 GeV at several values for the coupling strength; all curves are 
self-consistent final results. The solid lines show the result at the original coupling of G = 2.14 • A~^, the dashed lines 
were obtained for G = 2 • 2.14 • A^^ and the dotted lines display the width at a coupling strength of G = 4 • 2.14 • A^^. 
The larger couplings lead to a significant increase of the width. Since G enters the self-energies as a factor of G^ it is 
not surprising that a two times larger coupling results in a four times larger width and a four times larger coupling 
causes a width 16 times larger. The shape of Fq does not change very much. Only the gap in the antiquark sector 
gets smeared out due to the increased width which weakens the restrictions for the scattering processes. 

Figure [3 displays Fq for the different chemical potentials ojp — 0.268 GeV and ujp = 0.387 GeV at the regular 
coupling strength. The higher chemical potential has a similar effect as an increased coupling. Essentially, a scaling 
of the width can be observed while the general shape of Fq remains unaffected. For po > top one would expect a 
width Fo ~ p. This is, however, hard to see because of the vicinity to the cutoff in our calculations. For po < t^F, i-e., 
for the hole width, one expects a width Fq ~ (see Fig. (21 second line), because both momenta k and q lie within 
the Fermi sea. The observed scaling, however, is slightly smaller. This is due to the fact that the cutoff suppresses 
outgoing states with energies rp > A that can be reached when up > K/2. So the widths in our results differ not by a 
factor of 9 but only by a factor of 6 as one would expect for a Fermi energy oi uop = A/2 (this is the largest Fermi 
energy possible for which no outgoing states are suppressed). 

The spectral function defines the spectrum of possible energies po for a particle or antiparticle with momentum p 
that is added to the medium. In Fig. |S1 po is shown for the three coupling strengths at the lower chemical potential 
oi LOp = 0.268 GeV. The structure is clearly dominated by the two peaks at po — |p| and po = —\p\- These are the 
on-shell peaks of the quarks and the antiquarks. Because of the small widths they are very narrow and most of the 
strength is concentrated there. As one would expect from the results for Fq the peaks get broader and less high when 
the coupling is increased. Thus strength is removed from the peaks to the off-shell regions. The background of the 
spectral function seems to scale with the coupling in a similar way as Fq. Approximately a quadratic dependence is 
found. 
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The spectral function for the quarks and the antiquarks is not symmetric. This is due to the finite chemical 
potential: The quarks of our medium fill up all states up to the Fermi energy. Due to Pauli blocking only scattering 
processes with outgoing quarks at energies above cup are allowed. Incoming quarks from the medium must always have 
energies below LUp. On the other hand, all antiquark states are empty. Scattering processes with outgoing antiquarks 
at all energies are possible. However, there are no processes with incoming antiquarks from the medium. 

Since we are interested in short-range correlations the spectral function at large momenta above the Fermi momen- 
tum is of particular interest. The cut at \p\ = 0.3 GeV in Fig. |S1 indicates that short-range correlations seem rather 
weak for the original coupling strength, the spectral function is very small at < po < ^f- Only when the coupling 
is increased a significant population of high momentum states can be observed. 

In Fig. 1^1 the influence of the chemical potential on the spectral function is illustrated. As expected from the 
results for the width the peaks get broader and less high for the larger chemical potential of ojp — 0.387 GeV. Of 
particular interest is the cut for \p\ — 0.5 GeV (i.e., p > pp). A bump in the region of the occupied quark states 
indicates the growing importance of short-range correlations for higher ujp (one should compare this bump to the cuts 
for \p\ = 0.3 GeV in Fig. |Hlto take into account the different Fermi momenta). 

Finally we take a look at the (normalized) momentum distribution of the quarks which is given by 



The momentum distribution of nucleons in nuclear matter or nuclei shows a depletion of occupation probabilites by 
about 10% 5]. The resulting high- momentum tail is taken as a universal sign of short-range correlations. Figure 
1101 displavs the quark momentum distribution at the different coupling strengths for u)p = 0.268 GeV. As expected 
for an infinite system a sharp step at the Fermi momentum appears. At the lowest coupling a depletion of only 
0.1% is found. This confirms the previous interpretation of the spectral function. For the coupling twice as large 
the short-range correlations increase but still the depletion effect is below 1%. Only for the largest couplings we find 
a high-momentum tail of a few percent, comparable to the case of nucleons. Figure ITTl shows again that the larger 
chemical potential has an effect similar to the increased couplings. The depletion for ijjp = 0.387 GeV grows by almost 
one order of magnitude compared to the case of LOp ~ 0.268 GeV and the original coupling strength. It has a size of 
almost 1% and is comparable to the result for the two times larger coupling. 



In this paper we have presented an approach to the spectral function of quarks in quark matter. Based on a 
successful concept for nuclear matter it was shown how to construct a simple model for the quarks. Using the 
relations between collisional self-energies and the spectral function an iterative method was derived that goes beyond 
the quasiparticle approximation. Assuming that the exact structure of the interaction is irrelevant for the calculations 
- an important finding for nuclear matter - the pointlike interaction of the NJL model was used in the calculations. 

Our calculations, which are currently restricted to zero temperature and the chirally restored phase, show that we 
are able to deal with the numerical difficulties implied by this approach. However, the results also indicate that the 
infiuence of short-range correlations is small compared to nuclear matter. This finding, however, might be an artefact 
of the present model, the NJL model with vacuum parameters in the Born approximation. Since we know now that 
the model is technically feasible we can go beyond the simple model presented here. It is our plan to use more 
sophisticated interaction models for future calculations and to explore also other phases with broken symmetries. 
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FIG. 1: Lowest order contributions to the sclf-cncrgy iE^(p) in the SU(2) NJL model (Born diagrams). Since we have a 
pointlike interaction the direct and the exchange diagrams look the same in this figure. 
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FIG. 2: Scattering processes corresponding to the self-energies S< (p) in the Born approximation. The double lines carry the 
four-momentum p of ; p, q, k, r denote quarks with four-momentum {po,p), etc. (with po,qo,ko,ro > 0), and p, q, k, f denote 
antiquarks with four-momentum {—po,p), etc. (with po , qo , ko , ro < 0). Total collision rates arc found by integrating over the 
four-momenta k,q,r (respectively k,q,f), i.e., all lines except the double line. In the first row the scattering-out processes 
contributing to the loss rate lE^ are shown. The second row shows the only contribution to the gain rate — iE^ for positive 
Po- In the last row the processes contributing to the antiparticle loss rate — iE^ at po < are shown. Note that there are no 
contributions to E^ for po < ujf and to for po > u!f at zero temperature. 




FIG. 3: Examples for diagrams that are resummed in our self-consistent calculation. In this figure the lines correspond to 
free propagators. Once the elementary sunset diagram on the left is constructed from full propagators the other diagrams are 
automatically resummed. 
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FIG. 4: Series of diagrams that correspond to the coupling of dynamically generated mesons (a) and diquarks (b) to the quarks 
(see [Til. for details). The dashed line is a pion or sigma and the dotted line is a diquark. In the present model only the 
lowest order diagram, the sunset diagram, is included. 
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FIG. 5: Width of the spectral function of quarks at different momenta. The dashed Hne shows the result after the first iteration, 
the solid line shows the final result after two iterations. 
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FIG. 6: Width of the spectral function of quarks at different momenta. The solid line shows the width at the usual coupling 
strength {G = 2.14 ■ A~^), the dashed line is the result at a coupling twice as large and the dotted line has been obtained with 
a coupling four times larger than the usual value. 




FIG. 7: Width of the spectral function at two different momenta. The solid line shows the result at a chemical potential of 
ujf = 0.387 GeV while the dashed line has been found for ujp = 0.268 GeV (cf. Fig. [SJ. In both cases a coupling of G = 2.14-A"^ 
has been used. 




FIG. 8: The spectral function of quarks at different momenta and couplings (cf. Fig. El for details). 




FIG. 9: The spectral function for two different chemical potentials and momenta (see Fig. |7|for details). 
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FIG. 10: Quark momentum distribution in quark matter at different coupling strengths (see Fig. |S|for details). 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 



p [GeV] 

FIG. 11: Quark momentum distribution in quark matter for two different chemical potentials (see Fig. |3for details). 



